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Abstract 

We explain the relation between certain random tiling models 
and interacting particle systems belonging to the anisotropic KPZ 
(Kardar-Parisi-Zhang) universality class in 2 + 1-dimensions. The link 
between these two a priori disjoint sets of models is a consequence of 
the presence of shuffling algorithms that generate random tilings un¬ 
der consideration. To see the precise connection, we represent both 
a random tiling and the corresponding particle system through a set 
of non-intersecting lines, whose dynamics is induced by the shuffling 
algorithm or the particle dynamics. The resulting class of measures 
on line ensembles also fits into the framework of the Schur processes. 
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1 Introduction 


In this contribution we start by considering the domino tilings of the Aztec 
diamond, probably the most well-known and studied random tiling model, 
and we explain in quite some details how this model is equivalent to an 
interacting particle system. After that two generalizations are presented. 

The connection between the Aztec diamond and a (1 -I- l)-dimensional 
interacting particle system is well-known since the work of Jockusch, Propp 
and Shor [19], where they used it to prove the arctic circle theorem. Indeed, 
one quarter of the arctic circle is nothing else than the hydrodynamic limit 
of the height function associated with a totally asymmetric simple exclusion 
process in discrete time with parallel update and packed initial condition. 

In [7| we constructed a class of (2 -|- l)-dimensional interacting particle 
systems, some of which have the property that their fixed time projections 
can be viewed as a random tiling model, while certain (1 -|- l)-dimensional 
projections represent models in the Kardar-Parisi-Zhang universality class 
of stochastic growth models (like totally asymmetric simple exclusion pro¬ 
cesses). The domino tilings of the Aztec diamond is a special case that hts 
in that framework as remarked at the end of Section 2.6 of [7|. This was not 
a complete surprise, as the connection between the shuffling algorithm of the 
Aztec diamond and an interacting particle system in 2 -|- 1 dimensions was 
previously obtained by Nordenstam [23] (this was later extended to a few 
other cases in 121 ). 

The particle conhgurations in the state space of the (discrete time) 
Markov chain in [23] are not in a one-to-one mapping to domino tilings of 
the Aztec diamond, rather in a one-to-many relation. Here many is not an 
arbitrary quantity, but it is the number of random hllings in the shuffling algo¬ 
rithm during the last iteration. The mapping is the same as the one-to-many 
mapping between tilings of the Aztec diamond and the six-vertex model with 
domain wall boundary conditions at the free-fermion point [T8l|3T] . In order 
to have a bijection one needs to have the information of the particle conhg¬ 
urations at the two latest time moments (this fact is used for instance in the 
java applet m)- The reason for this fact will be clear when discussing the 
parallel update in the Markov chain, see Section |3l 

To explain the relation between the Aztec diamond and our interacting 
particle system, we map both systems to a set of non-intersecting (broken) 
lines viewed as paths on a Lindstrom-Gessel-Viennot (LGV) planar directed 
graph. They were hrst introduced by Johansson in [2U] in the study of the 
Airy 2 huctuations of the border of the disordered phase. The previously 
studied line ensembles representations as in [3T] are still one-to-many, while 
the one used in [T8l[2n] is in bijection with the tiling conhgurations. This 
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Figure 1: LGV graphs and Schur process. The weight of the left one is 
5A//i(«), the weight of the right one is 


bijection results in the shuffling algorithm [I61I29] inducing a Markovian dy¬ 
namics on the ensembles of lines. The same line ensembles with the same 
dynamics can be also obtained from the conhgurations of a system of inter¬ 
acting particles in 2 -|- 1 -dimensions in discrete time with parallel update [7] 
and with a particular initial condition. 

The probability measures (on tilings, sets of particles, or line ensembles) 
one obtains in this way fall into the class of so-called conditional L-ensembles 
and thus have determinantal correlation functions [9]. These measures and 
the corresponding Markov dynamics can also be described in terms of dimer 
models on Rail Yard graphs mm- 

Connection with the Schnr processes 

The constructions of the present paper can be recast into the framework 
of the Schur processes (introduced in [26]) and Schur dynamics on them 
(introduced in [I]). While we won’t do that here (see, however, a recent 
exposition in 0), let us make the connection more explicit. The basic fact 
is that the weights that we associate to the two sections of our LGV graphs 
with weights of vertical edges in the hrst graph being a and weights of skew 
edges in the second graph being jS (all other edges have weight 1 ), can be 
viewed as values of the skew Schur functions s\i^{a) and sa/^(/ 3). Here a 
stands for the specialization of the symmetric functions into a single non¬ 
zero variable equal to a (with complete homogeneous symmetric functions 
specializing into h„(a) = «"■, n > 0), and (3 stands for the specialization into 
a single ’dual’ variable equal to (3 (with complete homogeneous symmetric 
functions specializing into hn{j3) = 1 if n = 0, hn{/3) = f3 ii n = 1, and 
hnW) = 0 if n > 2). The commutativity relations of type fl3.18p below then 
turn into instances of the skew-Gauchy identity (see e.g. m)- 
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We chose to use more general notation in this paper (and in [7j) because 
of a potential applicability of the general setup to other models, which was 
indeed realized for example in [5]. 

Most of the work reported in this note was done in 2010, but it seems to 
not have lost all of its novelty yet. Indeed very recently we have been asked 
about the details of the connection between our work [7] and random tiling 
models. Thus we decided to make them available with this publication. 

Outlook 

In Section [2] we dehne the Aztec diamond, its domino tilings, and introduce 
the mapping to non-intersecting line ensembles. In Section [3] we present a 
class of Markov dynamics on interlacing particle systems and explain which 
initial conditions and parameters one has to choose to generate the same 
measure as the non-intersecting line ensembles of the Aztec diamond. Finally, 
in Section 0] we present two generalizations of the particle dynamics and 
explain which random tilings they correspond to. 
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2 The Aztec diamond 

In this section we introduce the Aztec diamond, a well-known random model, 
and explain his representation as a system of non-intersecting lines. The 
latter will be used later to make the connection with a Markov chain dynamics 
on interlacing particle conhgurations. 

2.1 The random tiling model 

Domino tilings of the Aztec diamond were introduced in [16]. For any A G N, 
we dehne the Aztec diamond A^r of size N as the union of all lattice squares 
[m,m -|- 1] X [n, n -1- 1], m,n G Z, such that they are inside the region 
{{x,y) : \x\ + li/l < A-|- 1}. A tiling of the Aztec diamond is a conhgura¬ 
tions of dominoes (i.e., 2x1 rectangles) such that the Aztec diamond 
Aat is fully covered by dominoes. Alternatively, one can think of a tiling as 
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Figure 2: Left panel: The Aztec diamond of size N = 10. Right panel; The 
associated bipartite graph. 


a perfect matching on the bipartite graph obtained by considering the dual 
graph inside the Aztec diamond, see Figure [2J Denote by Tn the set of all 
possible tilings of A^. It is well-known [16] that ITtvI = A random 

tiling of the Aztec diamond Aat is obtained dehned by choosing a probability 
measure on Tn- Specihcally we consider the probability measure given 
by 

I^SgTn ® 

where v{S) is the number of vertical dominoes in S, and a > 0 is a parameter. 

As one can see from Figure [3l a typical random tiling of the Aztec di¬ 
amond has four ordered regions and one disordered region in the middle, 
that has been studied in great details: the law of large number (Artie circle 
phenomenon) [T^], the fluctuations of the boundaries (Airy 2 process) [2U] . 
local statistics [13], and Gaussian Free Field fluctuations in the bullqj, and 
the GUE-minor process at the turning points 1211 . 

The shuffling algorithm 

To generate a random tiling of the Aztec diamond of size N, one can proceed 
iteratively from the Aztec diamond of size A^ — 1 by means of the shujfling 

^The detailed computations of the convergence to the GFF have not been written down 
in papers. The precise statement can be found in Section 6 of m- They are essentially 
the same computations as the ones made in Section 5 of |7], except that the kernel is the 
one associated with the discrete time TASEP with parallel update [S]. At the end of the 
argument, one uses ideas that goes back to m- A detailed statement of the expected 
result is available in Section 6 of m- 
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Figure 3: A random tiling of an Aztec diamond of size N = 200 and a = 1. 
The colors represents the four types of dominoes. 


algorithm introduced in [TB], that we recall briefly. For that purpose, let us 
give some names to the different types of dominoes: by superimposing the 
Aztec diamond to a checkerboard table (which of the chosen two possible 
ways is irrelevant), there are actually four types of dominoes, as the black 
site might be either to the right/left (resp. top/bottom) site of a horizontal 
(resp. vertical) domino. We call them North, East, South, and West domi¬ 
noes according to the following rule: the North (resp. South) domino is the 
horizontal one that hts into the top (resp. bottom) most part of the Aztec 
diamond, and similarly the East (resp. West) domino is the vertical one that 
hts into the right (resp. left) most part of the Aztec diamond. 

The shuffling algorithm is the following. Start with a domino tiling of A^v 
distributed according to Pat: 

Step 1: Move all the dominoes in the direction of their names by one unit. If in 
doing it a North and a South domino (resp. a East and a West domino) 
exchange their positions, remove them. 

Step 2: At this point, the dominoes partially tile Ajv+i. The empty region 
can be uniquely decomposed into 2x2 blocs. Independently of each 
other, each bloc is tiled with two horizontal dominoes with probability 
1/(1 -|- a^) or two vertical dominoes with probability a?/{I + a^). 

This procedure generates a random tiling of the Aztec diamond Aat+i with 
distribution Ptv+i as proven in [I6] . 
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Figure 4: A random tiling of an Aztec diamond of size A^ = 10 and its 
associates set of lines. The North dominoes are red, the South are yellow, 
the East are green, and finally the West dominoes are blue. 


2.2 Line ensembles for the Aztec diamond 

To each tiling of the Aztec diamond one can associate a set of non-intersecting 
lines bijectively as shown by Johansson in [20]. We use a slightly different but 
equivalent representation. Add horizontal lines in the middle of the South- 
type dominoes, a 45-degrees oriented segment in the West-type dominoes 
and a step-down in the East-type dominoes as indicated in Figure 01 For a 
tiling of Aat this results in N non-intersecting lines with fixed initial and end 
points. 

Further, one can think of the Aztec diamond of size N as being embedded 
into tilings of M^, where outside the Aztec diamond we add only horizontal 
dominoes that do not overlap and £11 the whole space, see the gray tiles in 
Figure [5] (left). By a simple geometric transformation, which can be easily 
recovered by comparing the left and the right illustrations in Figure [5l we 
obtain the final non-intersecting line ensemble representation of the random 
tiling of Figure 01 
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Figure 5: A random tiling of an Aztec diamond of size N = 10. Left panel: 
The original tiling of Figure 0] with the lines extended deterministically out¬ 
side A]s[. Right pane: Final representation of the ensemble of non-intersecting 
lines. The dotted lines is the underlying LGV graph. The lines which are 
hxed and correspond to the gray dominoes of the left picture are dashed. 


The above procedure gives a bijection between a tiling of an 
Aztec diamond and a conhguration of non-intersecting lines on 
a Lindstrom-Gessel-Viennot (LGV) directed graph. The basic 
building bloc of the LGV graph is the one here on the right. 

Then, if we consider an Aztec diamond of size N, we have a 
one-to-one bijection between the set of non-intersecting lines on 
the LGV graph that consists of N copies of the basic building 
bloc and where the lines start and end at all vertical positions 
— 1, —2,..., see Figure O for an example. 

Now we have a bijection between a set of non-intersecting line ensembles 
and the set of all possible tilings of the Aztec diamond A^r. The next question 
is to know which weights we have to assign to the edges of the LGV graph 
in order to obtain the distribution on A^r given by fl2.ip . i.e., where each 
vertical domino has weight a and each horizontal domino has weight 1. It is 
easy to see that it is enough to give weight 1 to each horizontal edge of the 
LGV graph and weight a to each of the vertical edges and the down-right 
edges. Alternatively, since the number of East and West dominoes are equal, 
we can also give weight 1 to each horizontal edge, weight a to each vertical 
edge and (3 to each down-right edge, provided that a[3 = a?. 
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Figure 6: Line ensembles decorated with domino colors as in Figure IH 

Induced line ensembles dynamics 

Because of the bijection described above, the shuffling algorithm induces a 
Markov dynamics on the set of non-intersecting line ensembles. To see how 
the dynamics works, it might be useful to decorate the line ensembles with 
dominoes of the same colors as in Figure IH see Figured The configuration 
at time n = 0 consists only of straight lines at height {—1, —2,..Given a 
line ensemble corresponding to an Aztec tiling of size u, the line ensemble of 
a random tiling of An+i induced by the shuffling algorithm is the following 
(see Figure [7] for an example): 

Step 1: Update the lines at the even lines-times (=horizontal coordinates) t: 
their height (=vertical coordinates) is updated to be equal to the height 
at lines-time t — 1. This is because by Step 1 of the shuffling algorithm, 
in the representation of Figure [6l the blue dominoes stay put, the red 
dominoes move vertically by 2, the green dominoes moves to the right 
by 1 and vertically by 1, the yellow dominoes moves to the right by 1, 
and when the trajectories of a red and a yellow (resp. green and blue) 
dominoes meet, they are deleted. 

Step 2: Update the lines at the odd lines-times t E {1, 3,..., 2?7, -|- 1}; if a line 
can be extended vertically by one unit without generating intersection, 
then this happens with probability /(1 -|-a^) (it corresponds to adding 
a blue/green domino, the vertical ones in the Aztec representation), 
while the move is not made with probability 1 /(1 -|-a^). All the possible 
updates in this step are independent of each other. 


9 






























o— 

o— 

r---Oc;-o— 

-- 

-o— 

----Oc;--:=0- 

r---Oc;-O 

----Oc;- 













Example of lines 

o—-^— 

for a random ti 

r---Oc;-O- 

----Oc: — — 

ing of Ai 

r - - - -Oc; -O 

----Oc;--^X) 


----Oc:-- 








• — 

- • -- 

- • -- 

- • - ^ 


After Step 1 



Figure 7: An example of the dynamics on the non-intersecting lines. The 
red dots are the points which, due to the shuffling dynamics, have to move 
(further indicated by a small arrow) or can more. After Step 1 the set of 
lines does not overlaps with the LGV graphs. This is the reason why the top 
line in Step 2 has to further move up. 
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Figure 8 : Graphical representation of GT 4 . The white dots represents the 
vector at level n = 4. 

3 Markov chains on interlacing particle 

In this section we hrst introduce a class of discrete time Markov chains on in¬ 
terlacing particle systems. It is a special case of the more general framework 
developed in [7]. With a particular choice of transition matrices and initial 
conditions, we can use the particle dynamics to generate non-intersecting line 
ensembles with the distributions equal to the ones induced by the shuffling 
algorithm. For pedagogical reasons, we focus here on the situation corre¬ 
sponding to the Aztec diamond. Two generalizations will be presented in 
Section 01 

3.1 Construction of the Markov chain 

The state space 

The state space of our Markov chain consists of interlacing conhgurations of 
particles, also known as Gelfand-Tsetlin patterns, dehned as follows: 

GTn = {X^ = {x\...,x^),x^ = K,...,x”) G ^ x^+\l<n<N-l} , 

(3.1) 

where 

-< ^ <x^< <...<<< (3.2) 

If x" -< we say that x" interlaces with We can (and will) think of 

the conhgurations as unions of levels: so the vector x” of is the state at 
level n. See Figure [ 8 ] for an illustration. 

The Markov chain is built up through two basic Markov chains: (a) the 
hrst one is the time evolution at a hxed level and (b) the second one is a 
Markov chain on GTat, linking level n with level n — l,2<n<N. We hrst 
discuss the two chains separately and then dehne the full Markov chain. 


11 


Markov chain at a fixed level 


The dynamics at a fixed level is a discrete time analogue of the famous 
Dyson’s Brownian Motion that appears in random matrices. Consider the 
one-particle transition matrix with entrie^ 


P(x.!/) = — 


. 1-p + pw ^ I _ 

dta--= < 1 — p, n y = X, 


To ^ " 10 , 

Then, the one-particle transition probability is given by 




ii y = X + 
iiy = x, 
otherwise. 


(3.3) 


Pt(x) = (P')(0,x) = CV'tl -p)‘ "llioj.Qi (3.4) 


for t G N. 

For any fixed n G {1,..., iV}, let us define an n-particle process on the 
projection of the state space GTat to the level-n, that is the Weyl chamber 

Wn, 

Wn = {X = (Xi, . . . , X„) G Z" I Xi < X 2 < . . . < Xn}. (3.5) 

Define also the boundary of the Weyl-chamber by 

dWn = {x = (xi,..., Xn) G Z"' I Xj = Xj+i for some i G {1,..., n — 1}}. 

(3.6) 

Consider first the process on x"' = (xj’',... ,xjj) G Z"' where each of the x^, 
1 < k < n, are independent processes with transition matrix P given above. 
The (discrete time) generator of the random walk is 

n 

(L„/)(x) = ^p(/(x-fei)-/(x)), (3.7) 

i=l 

with Cj the vector with entries ej(j) = We want to condition the random 
walk x" on never having a collision between any two of its components, i.e., 
to stay in the Weyl chamber forever. The conditioning can be achieved by 
the Doob h-transform with h given by the Vandermonde determinant, see 
e.g. 123. Indeed, one can verify that 

hn{,x) = (xj - Xi) = A„(x) (3.8) 

^For a set S', by ^ dwf{w) we mean the contour integral where the contour can 
be taken to be any anticlockwise oriented simple path containing all the points in the set 
S but no other singularities of the function /. 
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is harmonic and vanishes at dWn, i.e., L^hn = 0. Then, the conditioned 
process given that the walker never leaves the Weyl chamber Wn is the Doob 
h-transform of the free walk. More precisely, for x,y E Wn and t G N, the 
transition probability Pn,t from x to y oi the conditioned random walk is 
given by 

Pn,t{^, y) = ■Jj^^Prob(a:”(t) =y,T>t\ x"(0) = x), (3.9) 

where T = min{f > 0\x^{t) ^ Wn}- Since our random walks are one¬ 
sided, x^ can exit the Weyl-chamber only through dWn- Therefore using 
Karlin-McGregor’s type formula (or, since we are in discrete time, the LGV 
theorem), we have 

Prob(a:’^(t) =y,T > t\ a:''(0) = x) = det {pt{yi - Xj))P^-^. (3.10) 

Gonsequently, we define the level-n chain as follows. 

Definition 3.1. The transition probability of the Markov chain at level n is 
given by 

Pn{x, y) = det (P(xi, %))^=i (3.11) 

with P as in h3. lA) . x,y E Wn, and An the Vandermonde determinant. 


Markov chain at a fixed time 


Now we consider a Markov link between levels of GT^v that generates the 
uniform measure on GT^r given the value x^ of the level N, i.e., with 
x^ = (x)^ < X 2 < ... < x^) hxed. It is well-known, see e.g. Gorollary IA.41 
that 


ff of GTjv patterns with given x^ 


n 


N N 


j -i 




n 


Af-l 

n=l 


nl 


(3.12) 


Thus, the uniform measure on GTjv given x^^ can be obtained by setting 


„ , . Ar_i I AT. # of GT]y_i patterns with given x^ ^ 
Prob(x \x ) = —r— _ _^- rr —1 


ff of GTjv patterns with given x 


N 


Using fl3.12p we obtain 


(3.13) 


Prob(i'^ ‘ 11") = (IV - * 11| 




(3.14) 


Gonsequently, we dehne the following. 
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Definition 3.2. For any n G {2,...,N}, define the Markov link between 
level n and n — 1 by 


:= (n-1)! 




A,(x-) 


-Ir^n 


(3.15) 


where x" G Wn and x” ^ G WA-i- 

The uniform measure on GTjv given x^ can be expressed by 


■prV—1 

lln-i 


Af 


Ln=l 

AArfx 






(3.16) 


n=2 


There is an important representation of the interlacing through determi¬ 
nants. This is relevant when doing concrete computations. 

Lemma 3.3. Let x^ G and x^~^ G lLAr_i be ordered configurations in 
the Weyl chambers. Then, setting x^~^ = virt a “virtual variable”, we have 

= det(0(xf,xf“^))f^.^i, (3.17) 

with (j){x,y) = l[y>x] = ^ fro dw and 0(x, virt) = 1. 

Proof. The proof is quite easy. For x^“^ -< x^ one sees that the matrix on 
the r.h.s. of fl3.17p is triangular with 1 on the diagonal. Further, by violating 
the interlacing conditions, one gets two rows or columns which are equal, so 
that the determinant is equal to zero. □ 


3.1.1 The intertwining condition 

The key property used in our dynamics on GTm is the intertwining relation, 
illustrated in Figure [9l namely 

A^i := 2<n<N. (3.18) 

In our specific case, to see that fl3.18p holds one uses the Fourier represen¬ 
tation for Pn (see fl3.3p l and of A”_^ (see Lemma 13.3p . The intertwining 
relation can be obtained quite generically when the transition matrices are 
translation invariant, see Proposition 2.10 of [7] or Appendix |A] for more 
details. Intertwining Markov chains were first studied in H. 
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Figure 9: Intertwining of the Markov chains. 


3.1.2 The sequential update dynamics 

Now we dehne a discrete time Markov chain on GTthat we call sequential 
update, since the update occurs level by level and for the update of level n 
we use the updated values of level n — 1. In particular, the projection to 
the particles {x\,x\,...} is the sequential update dynamics of the totally 
asymmetric simple exclusion process (TASEP). 

Dehne the transition probabilities of a Markov chain on GTjv by (we use 
the notation = {x^{t),... ,x^{t))) 


P^{X^it),X^it + l)) = P,{x\t),x\t + l)) 

1 ^ Pk{x^{t),x^{t + l))h'l_^{x^{t + l),x^~^{t + 1)) 


(3.19) 


for X^{t),X^{t + 1) G GTat, see Figure fTOl for an illustration. 

One can think of P^ as follows. Starting from X^{t) = {x^{t),, x^{t)), 
we hrst choose a:^(f+l) according to the transition matrix Pi{x^ (t), x^ (t + 1 )), 
then choose x‘^{t + 1) using h+i)) ^ which is the con¬ 

ditional distribution of the middle point in the successive application of P 2 
and A^, provided that we start at x‘^(t) and hnish at x^(t + 1). After that 
we choose x^{t + 1) using the conditional distribution of the middle point in 
the successive application of P 3 and A 2 provided that we start at x^(t) and 
hnish at x‘^{t + 1 ), and so on. 

With our choice of P„’s and Ajj_^’s introduced above, the dynamics is the 
following: 


• x\ just performs a one-sided random walk (with jump probability p). 

• x\ performs a one-sided random walk but the jumps leading to x\ = x\ 
are suppressed (we say that x\ is blocked by xj). 
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Figure 10: Scheme for the sequential update dynamics. 


• X 2 performs a one-sided random walk but the jumps leading to X 2 = x\ 
are forced to happen (we say that x\ is pushed by x^). 

• Generally, performs a one-sided random walk (with jump probabil¬ 
ity p), except that it is blocked by x^“^ and is pushed by x'^z\ (whenever 
they exists). The updates are made from lower to upped levels. 


Conserved measures 

Given the Markov chain on GTjv described above, it is of interest to know 
which class of measures are conserved by the time evolution. 

Proposition 3.4 (See Proposition 2.5 of [7]). Let /iAr(x^) a probability mea¬ 
sure on ITat. Consider the evolution of the measure 

M^(X^) =/i^(x^)A)^_,(x^,x^-l)A))(:^(x^-^x^-2)•••A?(x^xl) (3.20) 

on GTtv under the Markov chain Pjf. Then the measure at time t is given 
by 





t times 


(3.21) 


For this statement to hold we use the intertwining property fid.isp . The 
proof is elementary. 
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An example of measure of the form of Proposition 13.41 is the packed initial 
conditions 

x^it = 0) = —n + /c — 1, 1 < k < n < N. (3.22) 

This initial condition is, due to the interlacing enforced by the A”_^’s, ob¬ 
tained by setting /iAr(x^) = It is interesting to note that such 

a measure has a determinantal form. 

Lemma 3.5 (See proof of Theorem 2.25 of [7]). Consider the probability 
measure on given by 

= const An{x^) (3.23) 

where 

Ti(x) = — / d^(l - l<i<N. (3.24) 

271-1 Jro 

Then kN{x^) = (5x^,(-Ar,...-p- 

This statement is a key for a more difficult Theorem 13.81 below. 

3.1.3 The parallel update dynamics 

Now we dehne a discrete time Markov chain that we call parallel update, since 
the update occurs in parallel for each level and for the update of level n we 
use the values of level n — 1 at the previous time instant. In particular, the 
projection to the particles {x\,X 2 , ...} is the parallel update dynamics of the 
totally asymmetric simple exclusion process (TASEP). 

The state space of the Markov chain is 

§JV = {x^ eWi,...,x'^ eWn\ Al_^{x^-\x^) >d,2<k<N}. (3.25) 

In our case, is as GTexcept that now the particles are weakly interlacing, 
i.e., the strict inequalities in fl3.2p have to be replaced by inequalities. There 
is an analogue of Lemma 13.31 (that we are not using here), see [8]. Dehne 
the transition probabilities of a Markov chain on Sjv by (we use the notation 
X^{t) = {x\t),...,x^{t))) 

P^{X^{t),X^{t + 1)) = P^{x\t),x\t + 1)) 

Pk{x^{t),x'^{t + l))A’l_^{x^{t + l),x’^-\t)) 

Ah- 

for X^{t),X^{t -I- 1) G Sat, see Figure [U] for an illustration. 
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X^{t) " X^{t + 1 ) 

Ia* 

1 A? 

xHt) - i’(t + i) 

Figure 11: Scheme for the parallel update dynamics. 


With our choice of P„’s and A”_^’s introduced above, the dynamics is 
similar to the one in sequential update: performs a one-sided random 

walk (with jump probability p), with the constraints that x'^ is blocked by 
x'l~^ and is pushed by x'^Zi- However, this time, the update of the particles 
is parallel, i.e., hrst one checks which particles could jump because they are 
not blocked and then these free particles independently jump to their right 
with probability p, pushing, when required by interlacing, other particles at 
higher levels. 

The analogue of Proposition 13.41 for the Markov chain on with parallel 
update is the following. 

Proposition 3.6 (See Proposition 2.5 of [7]). Let pp^{x^) a probability mea¬ 
sure on Hat. Consider the evolution of the measure 

M^(X^) = p^{x^)A%_,{x^, x^-^)A^-_l{x^-\x^-^) ■ • • A?(x^ 

(3.27) 

on Sat under the Markov chain . Then the measure at time t is given by 


{Mm ■ ■ ■ P^){X^) 



t times 


(3.28) 


Parallel update - extended state space 

In order to obtain the connection with random tiling of the Aztec diamond, 
we need to consider two small extensions. First of all, the transition prob¬ 
abilities can be taken time-dependent as long as the intertwining property 
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Pz{t-2) , 

x^{t) + 


Ai 


Ai 


y\t) 


P2{t — 2) X^(t) -^2(^—1) 

y\t + 1 ) ^ 

A? 


x^{t + 1 ) 


A? 


^ Pi(t-l) x\t) Pi(t) ^ 

+ 1 )-^ + 1 ) 

Figure 12 : Scheme for the parallel update dynamics with extended space 
time and time-dependent transition probabilities. The intertwining condition 
reads, A"_iP„_i(r) = Pn{T)K^_^. 


fld.lSp holds. We introduce the ’’time”-variable r, which denotes the hori¬ 
zontal position in the scheme as in Figure [T^l where as a reference point, 
r = 0, is set such that Pi(0) is the transition of a:^(0) to x^(l). Secondly, we 
extend the state space. Instead of considering the measure 03.271) as the one 
induced by the measure hn on hFAr through the action of consecutive actions 
of consider it as the induced measure by the consecutive actions 

of Pn-i- The middle points of the chain (A”_^ o P„_i)(x", is 

denoted by y"'. The choice of using the index n rather than n — 1 is slightly 
arbitrary, but the reason lies in the fact that it is natural to use the same 
index corresponding to the same “time”-r in the generalizations (see next 
section too). Thus we consider the extended state space 

= {X^ G §jv, e §Ar_i I e Wk, y’^ ^x\2<k< N}. (3.29) 

where = {x^,..., x^) and = (y ^,..., y^~^). The induced one-time 
transition of this extended state space can be described as following: 

• the y-particles are deterministically updated: y^~^^{t -|- 1) = x’^it), 

• the x-particles are updated according to the parallel update fl3.26p . 

see Figure [T 2 I Also in this setting there is a measure that is conserved. 

Proposition 3.7. Let fijsf{x^) a probability measure on hFjv. Consider the 
evolution of the measure at time t = 0 given by 

Mx^)A^_,{x^,y^)P^.,{T){y^,x^-^) 

X A^-_l{x^-\y^-^)PN-2{r + l){y^-\ x^-^) ■ ■ • 

xAl{x\y^)P,{T + N-2){y\x^) (3.30) 
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under the Markov chain with parallel update in the extended state space. 
Then, the measure at time t is given by 


(/iArPAr(r) • • ■ Pn{t + t - 1))(x^) y’^)PN-l{T + t){y^, 

X A^-_lix^-\y^-^)PN-2iT + t + x^-2)... 

X A?(x^ f)P,iT + t + N- 2){f, x^). (3.31) 

The proof of Proposition 13.71 is almost identical to the one of Proposi¬ 
tions [3]1] and [3T] (see Proposition 2.5 of [7]). 


3.2 Line ensembles associated with parallel update 

The map from a configuration in the extended state space to a set of non¬ 
intersecting line ensembles is simple. The LGV graph is the one of Figure [5l 
First let us describe where the particles are located and then we will explain 
how to add the lines joining them. The x-particles are at positions 

{(2fc - < j <k<N} 

and the ^/-particles are at positions 


{i2k,y’;+^),l<j<k<N-l}. 


At this point we can join the top particles at every time with a line following 
the LGV graph, then join the top remaining particles and so on. This results 
in a set of non-intersecting line ensembles as the solid lines in Figured] (right) 
(the deterministically added dashed lines are not yet present, but can be now 
added). 

The above discussion concerns the dynamics only and therefore it is clear 
that to obtain the Aztec diamond distributed according to the measure fl2.ip 
one has to start with a special initial condition and also to choose the tran¬ 
sition matrices P’n('r) in a r-dependent way such that particles at level n can 
move for the first time at time moment n. The initial condition is the packed 
initial condition described above, see Lemma 13.51 and the coefficients in the 
transition matrices P’n('r) are chosen to be p = a?/{I + a?) for r > 0 and 
p = 0 for r < 0, i.e., with P in Definition 13.11 given by 


1 

P(r)(x, y) = T— '^>0, 

A. (Jj X ~t~ ^ 

PiT){x,y) = r < 0. 


(3.32) 


The dynamics on the line ensembles is then exactly the one described at the 
end of Section [21 Indeed, first of all the initial conditions are the same. Then, 
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given the values after n steps, the parallel dynamics hrst update determinis¬ 
tically the values of the y-particles and then of the x-particles. These latter 
ores try to jump to their right with probability p, which is non-zero only 
for levels below n -|- 1 (compare Figure [12]). The jumps are allowed only if 
the arrival conhguration satishes interlacing, i.e., non-intersecting in terms 
of line ensembles. 


Theorem 3.8. The probability measure on the line ensembles (i.e., on the 
extended state space) corresponding to the Aztec diamond of size N is given 
by 


const X det[Tf (xf )]5=i det[0(xf, 2/f )]5=i det[/(|/f - xf 

X demxt\yf-Xf=\ det[/(l/f-^ - (3-33) 

• • • X det[0(x^, det[/(|/2 - 


where= ^^§.^^dz{l — z)^ with p = a'^/(l + a‘^), 

and f{x) = pSx^i -l- (1 — p)dxp, and (j) given in Lemma \373. 

This result is simply obtained by plugging in Theorem 13.71 the explicit ex¬ 
pressions of the packed initial condition (see Lemma [3.5p . the determinantal 
form of the of Lemma [3.31 and of Pn (see Dehnition 13.11 with 03.31) and 

Remark 3.9. It is well known that a measure of the type of Theoreml^ is 
determinantal as it fits in the general framework of conditional L-ensembles 
developed in If (see for example Theorem f.l of for a continuous case 
analogue). From this one can start computing a lot of observables and, in 
particular, obtain the law of large numbers / central limit theorem type results 
mentioned in the introduction. 


4 Generalizations 

In this section we present a generalization of the Aztec diamond and another 
random tiling model. Both have non-intersecting lines representations with 
Markovian dynamics induced by a shuffling algorithm, which can be linked 
with the intertwining Markov chains discussed in Section [3l 

There are other random tiling models with non-intersecting line represen¬ 
tations. One example is the double Aztec diamond studied in [Il[2]. They 
might ht in the general framework, but it is not so straightforward as the 
underlying LGV graph is not translation-invariant. In that case, the transi¬ 
tion probabilities are not translation-invariant and thus Proposition IA.5I can 
not be used to verify the intertwining property. 
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^0 , 1 /^ y’^ x^ 




/3i O 2 (32 Ct3 /S 3 «Af /Stv 

Figure 13: LGV graph for the Aztec diamond distributed according to fld.ip . 
The vertical edges have weights ai = a?q^^,a 2 = ... ,aN = o?q^, 

the diagonal edges have weights /9i = l/q^^~^, {32 = ... ,{3n = 1/q, 

and all other edges have weight 1. 


4.1 Aztec diamond with non-uniform weights 

In the original paper of Elkies, Kuperbert, Larsen and Propp [16], they con¬ 
sidered more general weights on random tilings of the Aztec diamond, namely 

PAr(T) oc (4.1) 

where v{T) is the number of vertical tiles and r{T) is what they called rank. 
The rank is the minimal number of elementary moves (switch of a pair of 
vertical dominoes into a pair of horizontal dominoes, or vice versa) that is 
needed to generate the tiling T from the tiling Tq consisting of exclusively 
horizontal dominoes. The value of r{T) can be read off the line ensem¬ 
bles. For the Aztec diamond of size N, only the top N lines, call them 
G(T),..., IjsiiT), are not deterministically set to be straight lines (see for 
instance Figure [5]). The line ensemble of the tiling Tq consists of all straight 
lines. Then, r(T) = Ylin=i — £„(To)). This model was further 

studied in [iQiiig. 

From the line ensembles of the Aztec diamond of size N, it is easy to 
verify that the weight Pjv can be obtained from the LGV graph made by N 
copies of the basic graph but with non-homogeneous weights. More precisely, 
the resulting LGV graph is represented in Figure [131 

This line ensemble is the one generated by the Markov dynamics with 
parallel update, with packed initial condition, and the transition matrices 
given as follows. Gonsider parameters ai, 0 : 2 ,... all different (the case where 
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(4.2) 


two or more are equal is obtained as a limit): for x”, g Wn, 


n+1 


with 

where 


... +1. ^ detjg,^ det[^(x7 - 

detlaf n”.i(l +/3»(^)/ttj) ’ 

^m,0 /3n(^)^m,— 1 ) 

/^n(r) = 


(3n for r > 0, 
0 , for r < 0. 

One can easily verify that, for r > 0, 


(4.3) 

(4.4) 


I 7/"+^—T 

= J Il]=i Pn ' , if -x] e {0,l},l<j < n, 
] otherwise, 


detiUx- 


(4.5) 

while for r < 0 is it just det[/„(x” - T)]i,j=i = ]Tj=i 

Further, the intertwining is satished for the following choice of (see 
Proposition IA.5P : for x" G IF^ and y'^ G PFn-i, dehne 


A^,(x^ 2 /") = 


det[a^^' det[fn{x^ - yj)]lj=i 

det[af n;=l(i - ’ 


(4.6) 


where yjj = virt is a virtual variable in order to write the formula in a more 
compact way, fn{x — virt) = a^, and 


Um) = J2<S., 


m.k' 


(4.7) 


k>0 


4.2 Random tiling of a Tower 

In this section we present a generalization that hts in our general Markov 
chain dynamics, then we will explain the induced dynamics on a set of non¬ 
intersecting line ensembles, and hnally we will explain how this generalizes 
to a shuffling-type algorithm for an associated random tiling model. 


Markov chains 

Now we consider a Markov chain on the scheme of Figure [TH It differs from 
the one of the Aztec diamond represented in Figure [T^bv the fact that the 
transition between x” and is a sequence of two A and one P transition. 
The middle points of the chain between x"' and are denoted by y"^ and 
respectively. Specihcally, let us consider the case where the transition are: 
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A(-l) 

a;^( 0 ) -" 

1 1 A‘ 

vHO) - 2/^(1) - 

1 Ai 1 Ai 

Pjf-i) 

z^(0) i"(0)=^(l) 

1 a; 

»'(•)) - 


A(0) 

^3(0) 

^3(0) 

Pi(0) 


0 :^( 1 ) 

2 /^ 1 ) 


Figure 14: Scheme of the Markov chain for the random tiling of a Tower. 


• from x^(t) to 

• from y'^{t) to kllz\, 

• from z'^it) to P 2 n- 2 it — n + l), 

where T’n('r) and A”_]^ is given as in fl4.2p and fl4.6p in the limit ai = a^, = 
... = a, 02 = 04 = ... = d, and jSi = [32 = ■ ■ ■ = f3. 

We consider packed initial conditions for the which by inter¬ 

lacing implies also 2 ;”( 0 ) = y"'{0) = a;"'(0), and dehne the evolution of the 
Markov chain is as follows: 

• update in parallel -|- 1) = 

• update y"'^^{t + 1) as the middle point of the two-step chain 
iP 2 n-i{t-n + 1) o Allzl)iy''+\t),x^{t)), 

• update + 1) as the middle point of the two-step chain 

{P 2 n{t - n + 1) O Al^_^){x^+^{t),y^{t + 1)). 

Theorem 4.1. The probability measure on the line ensembles (i.e., on the 
extended state space) corresponding to the Tower of size N is given by 


const X det[Tr(xf)]?A idet[0(xf,2/f)]J7=4det[0(yf,^f)]J^r’ 


X det[/(a;f-^ - det[cl>{xf-\ det[0(|/, 

N 


3 m,]-- 

N-\ V-l\]2V-3 


X • • • X det[/(a:,^ - ^J )]-^=1 det[0(a:,^ 2/j )]-^=1 JJ 2/”, (4-8) 


n=l 
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f2o 1/^ ^ 




a a (3 a a jS a a jS 

Figure 15: LGV graph for the a — a — (3 generalization. The vertical edges 
have weights a, a and the diagonal edges have weight [3. All other edges 
have weight 1. The lines start and end from positions flo = ~2,.. 


where ^ f ~ \^+/3/z), f(x) = (5x,o + /54,i, and 

(j) is given in Lemma l37S\ with = virt as well as Z2^_^ = virt. Further, 

2n 2n—l 2n—2 

i=l j=l k=l 

Non-intersecting line ensembles 

In a similar way as for the Aztec diamond, we can 
produce a set of non-intersecting line ensembles as fol¬ 
lows. The basic Lindstrom-Gessel-Viennot (LGV) di¬ 
rected graph of this generalization is the following, where 
the hrst column of vertical edges have weights a, the sec¬ 
ond column have weights a, and the diagonal edges have 
weights (3. Gonsider the graph obtained by N copies of 
this basic LGV graph as in Figure [T5l 

For a conhguration {y^,x^, z^,y‘^,x'^, ..., z^,y^,x^}, place the 
x-particles at positions 

{{3k — 1 , x^), 1 < k < N,1 < j < 2k}, 
the 2 /-particles at positions 

{{3k - 2,2/J), l<k<N,l<j <2k-l}, 
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Figure 16; Line ensembles decorated with domino colors as in Figure flSl 
below. 


and the z-particles at positions 

{{3k -3,z'^),2<k<N,l<j<2k- 2}. 

The ensemble of non-intersecting lines is the one where lines start (at line¬ 
time 0) and end (at line-time 3A^) from { — 1, —2,..The weight of such 
a conhguration generated by the dynamics described above is the following: 
vertical edges just before ^/-particles have weight a, the ones just after the 
2 /-particles have weight a, and the diagonal edges have weight (3. All other 
weights are set to be 1, see Figure [16] for an illustration. 

Random tiling model 

The set of non-intersecting line ensembles described above and generated 
by the Markov chain, can be associated to a random tiling model. For this 
generalization it is easier to think of a random tiling as a perfect matching of 
a bipartite graph, instead of a covering of a region. For the Aztec diamond, 
the N = 1 bipartite graph is illustrated in Figure [T7I (right). In the present 
case, the A^ = 1 basic building bloc is the one in Figure [T7| (left) consisting 
of a hexagon and two squares. As it looks like a small tower, we call this 
random tiling of a Tower. 

The random tiling model of size N consists in a perfect matching of the 
following graph: in the left-most column we have N basic blocs where the 
hexagons share a horizontal edge, this is followed by A^ -|- 1 basic blocs ob¬ 
tained by putting the left vertex of the hexagon at the right edge of the 
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Figure 17: The basic building bloc of the Tower (left) vs. the one of the 
Aztec diamond (right). One can think of random tiling of the plane with 
“dominoes” defined as two neighboring basic elements (squares, triangles, or 
kite for the Tower; only squares for Aztec) or, equivalently, think of perfect 
matching of the dual graph inside the Tower / Aztec. 


squares. This is repeated N times, so that in the end we have (3iV — l)iV/2 
hexagons. In Figure [18] we show an example for N = 5. The weight of a 
random tiling configuration is the following: each green domino has weight 
/9, each blue domino has either weight a or a depending its horizontal coor¬ 
dinate. 

Shuffling-type algorithm 

It is not too surprising that one can define a shuffling algorithm to generate 
the random tiling. For simplicity, consider the uniform tiling model. Start 
with a domino tiling of a Tower of size N. 

Step 1: Blue dominoes stay put. Red dominoes move up by one unit. Yellow 
dominoes move to the right and down by one unit. Green dominoes 
move right and up by one unit. Also, delete the dominoes whose tra¬ 
jectories intersect while doing these moves. 

Step 2: What remains is a partial tiling of the Tower of size iV -|- 1. The 
remainder can be uniquely decomposed into disjoint pieces, each of one 
is one of the following 5 configurations: (a) the basic building bloc, (b) 
hexagon plus a square (above or below), (c) hexagon, and (d) square. 
Then randomly tile all these pieces independently. 

The result is a tiling of a Tower of size iV-|-l that had the desired distribution. 
In Figure [19] we show a random tiling with uniform distribution of a Tower 
of size N = 100, corresponding to a = a = (3 = 1. 
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Figure 18: An example of random tiling of a Tower of size N = 5. There are 
two types of red/blue dominoes, depending on which a the corresponding 
column belongs to. 


By choosing the parameters a, a, and f3 not all equal to one, one can 
obtain different hgures. For instance, one can get a heart-like limit shape as 
in Figure [201 One can generalize the above setting to have m > 2 instead 
of 2 a-blocs. In that case, by choosing different a-values for the different 
weights, one can get for instance a limit shape with m — 1 fjords instead of 
one (see [T51I27] for similar behaviors in lozenge tiling). In the global scaling, 
one expects to see Gaussian Free Field fluctuations at the edges of the 
regular pieces one expects to see an Airy 2 process [201128] . and tip of the fjord 
one expects to see the Pearcey process 







Figure 19: An example of random tiling of a Tower of size N = 100 with 
parameters a = a = P = 1. For better visibility the background graph and 
the black dominoes are not shown. Also, the picture is rotated by 90 degrees 
with respect to the one in Figure [T8l 

















































Figure 20: An example of random tiling of a Tower of size N = 100 with 
parameters a = 1/4, d = 3/4, (3 = 1. For better visibility the background 
graph and the black dominoes are not shown. 
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A Toeplitz-like transition probabilities 

Here we give some results on transition probabilities which are translation 
invariant. They are taken from Section 2.3 of [7]. Our cases are obtained by 
limits when all Oj’s goes to the same value, say 1. 

Proposition A.l (Proposition 2.8 of [7]). Let Oi,..., a„ be non-zero complex 
numbers and let F{z) be an analytic function in an annulus A centered at 
the origin that contains all aj^ ’s. Assume that ^ 0 for all j. Then, 

for x" e Wn, 

det {ad y det(/(xr - yf))lj=i 


where 


F(ap)---F(a-^) 


= det (a,^T 


dz 


nz) 

^m+1 


(A.l) 


(A.2) 


A simple corollary for the specihc case of the transition probability in 
Dehnition 13.11 is the following. 

Corollary A.2. For F{z) = 1 — p -\- pz~^, it holds 


and 


a ^ ^ I A 

fim) = - (h dz -— 

27ri Jpo 


, it holds 


P, 

if m = —1, 


1-p, 

if m = 0, 

(A.3) 

0, 

otherwise, 


,n\\n 
j lli,j=l 

= A„(x"). 

(A.4) 


y’^eWn 

Proposition A.3 (Proposition 2.9 of [7]). Let oi, ..., a„ be non-zero complex 
numbers and let F{z) be an analytic function in an annulus A centered at the 
origin that contains all aj^ for j = 1,... ,n — 1 and assume that F{ay) ^ 0 
for all j = 1,..., n — 1. Let us set = virt and f{x — virt) = a^. Then 

71—1 

det {ad det{f{x^ - 2/r'))"i=i 


= det(a,^T (A.5) 


Fiay)---Fiay_,) 

A corollary concerning the transition kernel fl3.15p is the following. 
Corollary A.4. Let us choose F{z) = {1 — z)~^. Then 




>-m+l 


1, if m> 0, 
0, otherwise, 


(A.6) 
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and, with y'!^ ^ = virt, 


5 ; (n-l)!A„_,fe”-‘)det(/(i”-sr‘));j=i = A„(x”). (A.7) 

By the above results, we can define the transition kernels 


det(af)”. , F(a-,')-..F(a-:) 

\ I /2j = l 

for x"', e Wn, and 


(A. 8 ) 


1 1 

det (g/ det(/(x^ - 

detfajM"'. -^(®i 

\ I / 2,7 = 1 




(A.9) 


for x” G hhn and ^ G hhn-i- In our application, T„ is P„ and T^_i is 
The intertwining condition is then a consequence of the following result. 


Proposition A.5 (Proposition 2.10 of [7j). Let Fi and F 2 two functions 
holomorphic in an annulus centered at the origin and containing all the aj^ ’s 
that are nonzero at these points. Then, 


T„(Fi)r„(F2) = T„(F2 )T„(Pi) = rn(Pip2), 

T„{F,)T:_,(F,) = r”_,(F,)T„(f2) = r"_,(FiF,), 
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